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We propose that the entanglement of mixed states is characterised properly in terms of a proba- 
bility density function V(£). There is a need for such a measure since the prevalent measures (such 
as concurrence and negativity) are rough benchmarks, and not monotones of each other. Consid- 
ering the specific case of two qubit mixed states, we provide an explicit construction of V(£) and 
show that it is characterised by a set of parameters, of which concurrence is but one particular 
combination. V(£) is manifestly invariant under SU(2) x SU(2) transformations. It can, in fact, 
reconstruct the state up to local operations - with the specification of at most four additional pa- 
rameters. Finally the new measure resolves the controversy regarding the role of entanglement in 
quantum computation in NMR systems. 



Quantum entanglement is a unique resource for novel 
(nonclassical) applications such as quantum algorithms 
[f| , quantum cryptography , and more recently, metrol- 
ogy [3|. Thus, it has a pivotal role in quantum informa- 
tion theory. It is also central to the study of the foun- 
dations of quantum mechanics However, while pure 
state entanglement is well defined, mixed state entan- 
glement (MSE) is still rather poorly understood. Cur- 
rently used definitions such as entanglement of formation 
(EOF) @ and separability [6| are based on the empha- 
sis given to a particular quantum feature. These defi- 
nitions are not equivalent 0] and are operational in a 
limited sense. Thus, concurrence as a characteristic of 
EOF Q is defined only for a two qubit system; nega- 
tivity as a criterion for separability 0, is necessary 
and sufficient only for two qubit systems and a qubit- 
qutrit system. Likewise, majorization (TlT ] is a necessary 
condition for separability. Further, concurrence and neg- 
ativity are not relative monotones, although the former 
bounds the latter from above. In particular, states with 
the same negativity may have different concurrence and 
vice versa. Note that real systems are almost always in 
a mixed state. Indeed, NMR quantum computers (NMR 
QC) [H, [l3j are prepared in the so called pseudo pure 
states which are highly mixed. Their concurrence (and 
hence negativity) is zero, and yet nontrivial nonclassical 
gate operations with up to eight qubits have been re- 
ported [l3|. More recently, a 12-qubit pseudopure state 
has been reported for a weakly coupled NMR system [l5| . 
To unravel the sense in which the entanglement is a re- 
source in these systems, there is a clear need to go beyond 
the above mentioned benchmarks. We address this prob- 
lem here, and propose an alternative definition of MSE. 

To motivate our approach, we recall that a mixed state 
is required to describe an ensemble of quantum systems 
each of which is in a pure state. Entanglement has a 
sharp value in each pure state; Thus, MSE may be ex- 
pected to acquire a statistical character, and be char- 
acterized by a suitably defined probability density func- 
tion (PDF). We propose below a definition of MSE, in 



terms of one such PDF, which is strictly operational and 
applicable to any bipartite system. The definition does 
not require any new notion of entanglement other than 
that for pure states. We proceed to give an explicit con- 
struction of the PDF for the important case of two qubit 
systems. For these sy terns, wc show that the PDF has 
some striking morphological features which completely 
encode the information on MSE: these features appear 
as a few points of discontinuity of various orders in the 
PDF. These points are shown to allow an almost com- 
plete reconstruction of the state, up to local operations 
(LO). It is shown how concurrence gets reinterpreted as 
a benchmark. Finally, the issue of entanglement in NMR 
QC gets naturally resolved. 

We now posit a probability density function for entan- 
glement, V p (£). The definition will be given in several 
steps: Let the state p of a two qubit system be charac- 
terised by its eigenvalues A|, with respective eigenstates 
\ipi) which are orthonormal; the notation implies that 
the eigenvalues are arranged in a non increasing order. 
The choice of \ipi) is non unique if the eigenvalues are 
degenerate, but it is of no concern to us here, (i) As the 
first step, we define a sequence of projection operators 

i 

n * = 5^1^^' 1 5 Ul c ni+1 ' with Ui 

being the full 

Hilbcrt space. It is a trivial identity that 

P = (Ai - A 2 )n x + (A 2 - A 3 )n 2 + 

(A 3 -A4)n3 + A 4 n4. (1) 

The above equation resolves p into an incoherent sum of a 
hierarchy of the subspaces IL , with the weights given by 
the nonnegative vector A = (Ai— A 2 , A 2 — A3, A3— A4, A4), 
whose norm is a measure of the purity of the state. In a 
sense, the vector represents the manner in which the state 
"spills over" to successively higher dimensional spaces, 
(ii) As the next step, observe that if p is a projection IL of 
dimension i, the ensemble would be uniformly distributed 
over states in IL: (ip\p\ip) = lV|-y>) £ II;. A Probability 



2 



Density Function (PDF) may be naturally defined thus: 
/ / dS'dHi 6(8' - 8) 



J dHi 



(2) 



with dHi being the appropriate Haar measure, (hi) As 
the last step, rescale p — ► p s and rewrite it in terms of 
the difference in relative weights pi = (A, — Aj +1 )/Ai as 
p s = pj\i. The definition of the PDF is then given 
by the simple superposition 



V p {£) = Y,»iPi{£). 



(3) 



The rest of the paper is devoted to an elucidation of Eqns. 
2 and 3. We choose the pure state concurrence 2|aj j ay — 
at,i.aj.i-|, in terms of the coefficients of expansion of \ip) = 
a TT | ft) + «nl U) + a ll\ it) + au\ II), as the measure 
of pure state entanglement. Although the definition is 
given for the simplest case, the generalization to higher 
spin systems is straightforward, and we do not discuss it 
any further in this paper. 



TWO QUBIT PROBABILITY DENSITY 
FUNCTIONS: DESCRIPTION OF SUBSPACES 



V(£) 




FIG. 1: Some Typical probability density functions for II2. 
Note the solid curve, which shows all the features of Vi(£}. 
It has a cusp at £ CUS p = 0.8 and goes to zero at £ m ax = 0.89. 
The step function is an extreme example, where £ CU sp = 0, 



and the other dotted curve, has £ cus 
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(£ ma x, V 2 {£ ma x))- or equivalently, (£ cusp , V 2 (8 max )) for 
characterizing the curve. A straightforward computation 
establishes the relations 



We first consider the situation when p is a projection, 
case by case. We then move on to discuss the general case 
(displayed in Eqn.3). Since the normalization is provided 
by dividing by the total volume of the group space, the 
trace factors will be dropped. We employ LO on the 
subspaces freely, since the PDF remains unaffected. 
The pure state: Consider p — fix = \>fi)(4>\- The proba- 
bility density function V\(8) is simply 5(8 — 8$), in terms 
of the the entanglement of \(j>). The PDF is singular, and 
specified by a single number. 

Two dimensional projection: p = II2 is the most 
complicated and the most interesting case. Suppose 
IV') £ II2. Let \xi): IX2) be orthonormal and span U 2 . 
We have, \ip) = |xi)cos§e^/ 2 + \ X2 ) sin § e'^l 2 . The 
Haar measure is simply read off as dTL = sin 8d9d<fr. By 
a suitable LO, we can choose \xi) to be separable, in 
its canonical form (1,0,0,0) in a separable basis, i.e., 
Ixi) = I IT)- IX2) can be further chosen to be of the form 
(0, x, y, +yl — x 2 — y 2 ), where x,y > 0. The entan- 
glement distribution is, therefore, characterized by two 
non-negative parameters, and is implicitly determined by 
Eqn.2. 

The generic form of the PDF in II2 is shown in FIG. 
1 (the solid curve). We observe that it has three mark- 
ers, (i) £ C usp, the entanglement at which the probability 
density diverges, invariably as a cusp, (ii)8 max , the maxi- 
mum entanglement allowed, and (iii) V 2 (8 max ), the prob- 
ability density at 8 max - In fact, any two of them suffice to 
characterise the PDF completely. One may specify e.g., 



8 m ax = xy + \J z 2 + x 2 y 2 (4) 
z 2 

£cusp ~~p, £"inax COS /X (5) 

&max 

fl = S[li l \F vlF W 

■ -i f 2y/xy(xy8 max + z 2 )\ 
= sin Jyz ( fi ) 

\ <^max J 

which allow us to determine the parameters x, y that de- 
fine II2. p is well defined by virtue of the inequality, 
V 2 (8 ma x) > ^/8 ma x- Note that unlike with the other 
measures, the state itself can be reconstructed up to LO. 

Two extreme cases occur when 8 cusp — and 8 cusp = 
8 ma x- In the first case, the PDF is a step function, ter- 
minating at some 8 max . In the second case, the density 
increases monotonically, diverging at 8 max (FIG. 1). The 
relative abundance of entangled states is more in the lat- 
ter case. One may per se expect that the associated 
concurrence should also be larger. Interestingly, how- 
ever, the concurrence is related to the new parameters 
by C = {8 max - 8 cusp )/2, vanishing when 8 cusp = 8 max . 
In other words, it is not sensitive to the relative abun- 
dance at zero (or small entanglements) at all. In any 
case, C emerges as a particular benchmark of the prob- 
ability density, characterizing it only partially. We note 
that if p — II3 or II4, then its concurrence is zero. By 
the convexity of the concurrence, we conclude that the 
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concurrence C p is bounded by 

C P < (Ai-A 2 )C ni + (A 2 -A3)Cn 2 . 

Incidentally, the entanglement distribution of a subspace 
IIS; orthogonal to II2 is the same as that of IT2. 
Three dimensional projection: We now move on to 
the case p = II3, whose PDF has a simpler structure. II3 
is completely characterised by its dual, |£) _L IT3. Thus, 
the PDF is characterised by a single parameter £j_ , which 
is the entanglement of the orthogonal state 
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FIG. 2: A Typical probability density for II3. Note the point 
of discontinuity in the derivative at £ = £± 

The integrating measure may be conveniently writ- 
ten as dH.3 = sin 2/3 sin 20 sin 2 6dad/3djd0, when the 
state is expanded in an orthonormal basis as: = 
cos0|xi) +e l< - a+ ^ sin6>cos/?|x2) - e l( - a -^ sin6>sin/3|x3), 
with the integration ranges, 6,0 g [0, f ] and a, 7 € [0, ir}. 
Conveniently, one may choose Xi,2 to be separable, and 
by a suitable LO, they can be brought to the form 
I Tt)i I ID- We have verified that the resulting proba- 
bility density can be cast into the simple form 

of 1 
Vs(£)= cosh- 1 ! — ). (7) 

V 1 - £± b > 

where £ > — max{£,£±). 

A typical curve for V^{£) is shown in FIG. 2, which 
exhibits the required characteristic. The curve possesses 
a discontinuity in its derivative at £± . Significantly, con- 
currence (being identically zero) fails to distinguish dif- 
ferent three dimensional projections, e.g., £± = or 1, 
although their PDFs are vastly different. 

Lastly, we consider the full space II4, whose PDF is 
universal. This curve is obtained by using the Haar mea- 
sure on Note that the curve is smooth every- 
where, as shown in FIG. 4. 

It remains to consider the case when p is an incoherent 
sum of the projections (see Eqn.3), where the weights 



V(£) 
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FIG. 3: The probability density V±{£ ) for the entire Hilbert 
space. 

have been chosen such that the special cases p = II are 
naturally recovered; they also ensure that the results are 
not artefacts of any basis. If \\p\ — P2W is small, the 
corresponding probability density functions will also be 
close to each other. 
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FIG. 4: The overall probability density Va^S) for a typical 
mixed state, p, with eigenvalues {0.385, 0.288, 0.231, 0.096}. 
Note that the features of the indivdual subspaces are vividly 
preserved. 

FIG. 4 illustrates the PDF for this general case. The 
important point to be noted is that the superposition 
of curves does not obliterate the information contained 
in individual curves; they are retained as points of dis- 
continuity or singularity and each individual PDF may 
be reconstructed, together with the associated weights. 
V{£) is by definition invariant under LO. With this, one 
may ask if the state itself may be reconstructed, up to 
LO. Before we take up this question, we consider an im- 
portant application of this prescription, to NMR QC. 
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NMR QC employs the so called pseudopure states for 
computation. Since it is experimentally demonstrated 
that the quantum logic operations used in QC are imple- 
mentable with NMR, it follows that these states should 
possess a non vanishing entanglement. Indeed, they 
have the form p pps = ^j^l + e\ip)(tp\, in our system of 
expansion. The NMR signal is sensitive only to the pure 
component, the so called deviation matrix. Accordingly, 
its V p {£) is given by a weighted Dirac Delta superposed 
on the PDF coming from the full space. The uniform 
background is invariant under unitary operations, but 
the one dimensional fluctuation is not, allowing for 
non-trivial gate operations. Thus NMR QC exploits 
the excess of entangled states over the unpolarized 
background as a resource, and this feature is correctly 
captured by the PDF of the state. This is in contrast to 
other measures which attribute a zero entanglement to 
all PPS with e < |, while usually in practice in NMR QC 
e ~ 10~ 6 . This analysis also raises the interesting pos- 
sibility of QC with more general pscudo projection states. 

Lastly, we return to the issue of the reconstructibility 
of the state (up to LO). If p is a projection, the recon- 
structibility is assured, by construction. When p is more 
general, the reconstruction is partial. For, the action 
of SU(2) x SU(2) on p produces an orbit of dimension 
six, characterised by nine invariants. The set of param- 
eters which characterize the entanglement are seven in 
number (for example: {/^., fi 2 , M3, £b £cus P , £max, )• 
Geometrically, V p (£) is invariant under independent LO, 
Li, acting on the subspaces IF, where IF C IF + i. If p is 
to be unique up to a global LO, one needs the additional 



constraint Li = UUp , where Lf may be chosen freely. 
Let us choose L^ = 1 (where 1 is the identity opera- 
tor). The ncstcdness condition, viz., that G IF and 



(0) 



\ip4,) G II§ , entails that L^ and L^ 1 get specified by two 
parameters each fl8j | . 

More explicitly, if we have IF in the canonical 
form, it is spanned by and \\2) given respec- 

tively as: (1,0,0,0) and (0,x,y,z). Therefore, we 
can specify = cos fe 1 ^ 2 + \ X2 ) sin fe"^/ 2 

by giving the values of (0,<p). Similarly, IV'J.) 
can be specified by (0±,4>±) when it is expanded 
in the canonical basis of IIj = (1 — IF), given 
by \xl) = (0,0, c/Vc 2 + & 2 ,-VVc 2 +6 2 ) and 

\X°2 



(0) 



i) = (0, Vc 2 + b 2 , ab/y/c 2 + b 2 , ac/y/c 2 + b 2 ). 



In conclusion, we have given a prescription that de- 
scribes the entanglement of mixed states by not just a 
number, but an exhaustive set of parameters which char- 
acterize the manner in which the entanglcnment is dis- 
tributed over the ensemble. They further permit an al- 
most complete reconstruction of the state up to LO. The 
prescription may provide a better insight into other mea- 



sures of entanglement such as entanglement of distillation 
and entanglement cost. Investigations along these lines, 
and a further study of the PDFs for higher spins may 
provide us with a better appreciation of quantum entan- 
glement. 
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In fact we need only one parameter each to fix the states 
up to discrete ambiguities. This is because from the Vi{S) 
we already know the entanglement of the states, which 
fixes one parameter. However, there remains a discrete 
ambiguity if only one of 9 or <f> is specified 



